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In a recent paper we considered the possibility of a scalar field providing an explanation for the 
cosmic acceleration. Our model had the interesting properties of attractor-like behavior and having 
its parameters of O(l) in Planck units. Here we discuss the effect of the field on large scale structure 
and CMB anisotropies. We show how some versions of our model inspired by "brane" physics have 
novel features due to the fact that the scalar field has a significant role over a wider range of redshifts 
than for typical "dark energy" models. One of these features is the additional suppression of the 
formation of large scale structure, as compared with cosmological constant models. In light of the 
new pressures being placed on cosmological parameters (in particular Ho) by CMB data, this added 
suppression allows our "brane" models to give excellent fits to both CMB and large scale structure 
' data. 
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I. INTRODUCTION 



Current evidence that the expansion of the Universe is accelerating jlj], if confirmed, requires dramatic changes 
in the field of theoretical cosmology. Until recently, there was strong prejudice against the idea that the Universe 
could be accelerating. There simply is no compelling theoretical framework that could accommodate an accelerating 
Q\ , universe. Since the case for an accelerating universe continues to build (see for example attempts have been made 
t-H ■ to improve the theoretical situation, with some modest success. Still, major "fine tuning" problems remain. 

All attempts to account for acceleration [3-31] introduce a new type of matter (the "dark energy" or "quintessence") 
with an equation of state = w^p^ relating pressure and energy density. Values of w$ < —0.6 today are preferred 
by the data p2| ] and in many models can vary over time. (In this framework, a plain cosmological constant is the 
case where = —1 independent of time.) Most models with a varying are based on a simple scalar field with a 
specific potential in the Lagrangian. In addition to the standard ingredients fed into the Friedman-Robertson- Walker 
(FRW) cosmology, one introduces a scalar field which in general has a varying w^. If the value of 1 + 3W^il^, becomes 



negative (with the assumption that the photon and neutrino contributions to the total density of the universe is small 



compared to the rest of the mass-energy) the universe enters an era of accelerated expansion. 



The evolution of spatially homogeneous scalar fields in an FRW cosmology has a long history in the context of 
cosmic inflation, but the inflaton fields do not play a significant role today. In the last two years scalar fields have 
been considered which can produce an accelerated expansion in the present epoch This area has been stimulated by 
the growing evidence for cosmic acceleration today. The additional scalar field matter is known as "quintessence" or 
"dark energy" . In the early models one had to fine-tune the initial conditions of the field in order to get an accelerated 
expanding universe today, a feature which is very undesirable. Later on the very interesting category of "tracker" 
quintessence models were created, in which the field gives the desired behavior independently of initial conditions. 
Still one had to introduce a small scale into the Lagrangian in order to achieve this. One way forward with these 
models is to try and construct a specific explanation for these small parameters J2^]2!| . 

In a recent paper |[9[ we discussed a class of quintessence models which like the "tracker" models, would produce 
the desired effect of an accelerated universe, independently of the initial conditions. The model was based on the 



pure exponential potential which was known to possess attractor solutions 1 33 39 1 . However, the pure exponential 
however, in the attractor regime, cannot produce a realistic accelerated expanding universe. In our previous paper we 
showed that when the exponential potential is modified by a polynomial prefactor, we can keep the attractor solutions 
during most of the history of the universe (and therefore the independence from the initial conditions) while at the 
same time producing an accelerating scale factor today. The new nice feature of the model, not found previously in 
other models, was that all the parameters involved were 0(1) in Planck units. Since then, a number of authors taken 
this idea in interesting directions fif],|d| , which adds to the case for pursuing the observable effects of our model on 
cosmic structure. (We also note that yet another class models using simple Planck-scale physics have been proposed 

In this paper we consider the evolution of small perturbations for the former potential and discuss the implications 
on the Cosmic Microwave Background (CMB) and structure formation. In particular the recent Boomerang and 
Maxima data suggest that the Hubble constant might be larger than the so far accepted value of 65kms~ 1 Mpc~ l . 
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In this case if one normalizes the data to COBE, the standard A-CDM model doesn't do that well for the observed 
matter power spectrum. The former potential seems to do slightly better. One of our important results comes from a 
modification to the potential reported in [jl9| , inspired by "brane" physics which results in a number of novel features. 
One of these features is a much better joint fit to the matter power spectrum and CMB anisotropy power. 

The structure of this paper is as follows. First we review the evolution of scalar fields in an FRW cosmology. We 
discuss in particular the two mentioned potentials (the one from (l^] and the "brane" variation) emphasizing the 
relevant physics and noting the differences between the two of them and between other models. Then we present 
the evolution of small perturbations for the two fields at hand and discuss the various effects. The evolution of 
perturbations is different in the two models which leads to considerable differences in the power spectra. The effect of 
the additional dark energy in the universe on the matter power spectrum is discussed next. We consider parameters 
of the models that give good erg's and fit the newly published decorrelated data of Hamilton et.al. p2] , p3[ . Following 
that we focus on the related physics that affect the CMB and also discuss the difference between the two models and 
A models. A final section gives our conclusions. 



II. EVOLUTION OF THE BACKGROUND 



Let us review the background evolution for a scalar field cosmology. This section does not introduce anything new 
(except the final part which introduces the new potential) and is merely intended for the new readers not familiar 

with the subject. We use a (H ) metric signature, with Greek indices running from to 3 and Latin indices 

(spatial) from 1 to 3. Index free 3-vectors are denoted as x and 4-vectors as x. Unless otherwise told, all units will be 
Planckian with M p = -7= = (SnG)^ 1 ^ 2 ^ = 1 , k e being Einstein' constant and M p being the Planck mass. Moreover 
we assume that the zeroth order cosmology is flat FRW in the conformal synchronous gauge. 



A. Cosmological scalar field evolution 



The point of departure is the action functional for the scalar field and gravity, 

S[9ap, 4>]= J dW=ff{f - + V(4>)} (1) 

Since the background metric has been fixed by our assumptions we can use ^/—g = a 4 and R — —6-%. Moreover 
the assumption of homogeneity and anisotropy forces the spatial derivatives of the field to zero, <j> ^ = 0. The action 
simplifies, and takes the form of one particle Lagrangian mechanics in two dimensions with coordinates a(r) and </>(t). 

S[a(r), 0(r)] = J dr{3d 2 - l -a 2 ^ 2 + aV(0)} (2) 
The canonical momenta are p a = 6a and — a 2 4>. The constraint H — gives the first Friedman equation 

and therefore the field density is 

P = ^4? + V(<t>) (4) 
Variation w.r.t. the field <fi{ T ) gives the field equation of motion 

4> + 2-<j) + a 2 V 4> = (5) 
a 

where denotes derivative with respect to <p. 

Finally variation w.r.t. the scale factor a(r) gives the second Friedman equation 

6-^r = -^U 2 +W = p-3P (6) 
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from which we get the field pressure as 

p = ^ 2 - Vtf) (7) 

If we want our scalar field to accelerate the universe we need the field to provide the dominant form of energy 
density and at the same time have negative pressure. In order to accelerate the universe a necessary and sufficient 
condition is that the deceleration parameter q given by 

q= -™ = -(H-3to^ + fi r ) (8) 
be negative or equivalently P$ < —H 2 (In the last equation the dot denotes a derivative with respect to real time). 



B. The pure exponential potential 

One very interesting potential is the exponential V — Vbe~ A( ^. This potential has been shown to have attractor 
solutions J33|]3~i| ] , that is regardless of initial conditions, the field eventually scales like the dominant matter component. 
For a detailed discussion see ||^] . 

If the dominant component scales as p n — Po(^f)™ then the scalar field approaches an attractor solution and its 
fractional density is given by = w ■ This is a special case of scaling a field - a classification of such fields can be 
found in [?]. 
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FIG. 1. LEFT: Fractional densities as a function a for the exponential potential. The solid curve corresponds to matter 
with £l m ps 0.88 today, the dotted curve corresponds to radiation T = 2.726K today and the dashed curve corresponds to the 
scalar field with f2</, ~ 0.12(A = 5) today. Radiation-to-matter transition is around a — and a = 1 today. 
RIGHT: Evolution of the equation of state parameter w$ and the deceleration parameter as a function of a for the same 
model. The solid curve corresponds to w and the dotted one to q. The attractor causes to simply mimics the w of the 
dominant form of matter so q^ can never be negative in the attractor regime of a pure exponential potential. 



During the attractor era the relevant physical quantities are 

n na' 

■ = ( ■-' 
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n - 2 r 



(9) 
(10) 
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where n = 3(1 + w). 

Figure |l| shows the evolution of such a field with the scale factor. Notice the change of the fractional density in the 
field from the radiation to the matter era. Since w simply mimics the w of the dominant matter, it is impossible to 
get an accelerating universe in the attractor era. It is possible to avoid the attractor if one chooses A < y/n. In that 

i 2 

case the attractor is not there but the field will have its own scaling behavior with = 4j — 1. It is also possible 
to choose initial conditions sufficiently far from the attractor solution that the attractor has not been reached even 
today @. 

C. Modifying the exponential potential 

The attractor behavior of the exponential is appealing, and it would be nice if we could modify it in such a way as 
to keep this behavior for most of the history of the universe and at the same time get an accelerated universe today. 
Indeed we showed in a previous paper |l9| ] that this is possible by including a prefactor in front of the exponential 
V = V^(0)e~ A( ^. The effect of the prefactor is to introduce a local minimum in the exponential such that the field 
gets trapped into it. After the field gets trapped it starts behaving like a cosmological constant and the universe 
eventually enters an era of accelerated expansion. In our previous paper we used a polynomial prefactor. The full 
potential then takes the form 

V = V [(4>- B) 2 +A]e' x ^. (11) 

Through this paper we refer to this potential as the "AS" model. 

The left panel of Fig. || shows the cosmological evolution of densities in an AS model. One noticeable feature of this 
model is the breakaway from the attractor behavior at late times. This will be important for structure formation as 
the field density becomes very small when structure begins to form and so it doesn't affect structure formation much. 
Lastly we should note that trapping takes place independently of initial conditions of the quintessence field in the 
very early Universe. This is because the Universe is drawn into the attractor at early times. Since the field is in the 
attractor solution by the time it approaches the local minimum there is no memory of the initial conditions. There 
is a variety of possible behaviors at the local minimum. If the minimum is very shallow however the field can roll 
through rather than be trapped. It is also possible to choose parameters which remove the local minimum, but the 
field still lingers long enough to cause an era of accelerated expansion. For trapping a necessary condition is AX 2 < 1. 
The various behaviors and their corresponding parameter ranges have been investigated in [|44| . 

Instead of choosing a polynomial prefactor it is possible to have the desired accelerated expansion by using a 
different function V{(j>). An oscillating term for example could do the job pl| . Here we give another form which can 
accelerate the universeFI. 



V 



° D 



B) 2 +A 



e -A * (12) 



Such a potential could arise from the various brane models as a Yukawa-like interaction between branes p^ j (we will 
call it throughout this paper the "Brane" model). 

The right panel of Fig. |^ shows the evolution of densities with the Brane potential. During radiation era it behaves 
like the AS potential. During the matter era however it retains significant density unlike the AS potential. This has 
a significant impact on both structure formation (which is suppressed) and the CMB as we shall see later. 

To understand why there is such a difference we plot the equation of state parameter w with the scale factor in Fig. 
| In the case of the AS model, soon after the field enters the matter era attractor where w should go eventually to 
0, the field instead rushes toward the minimum of the potential and then undergoes damped oscillation. During that 
time the energy density of the field falls of significantly faster than the than energy density of the pressureless matter. 
This behavior corresponds with a brief rise in w(a). Due to damping from the expansion however, the oscillations 
eventually stop, the field settles in the minimum and starts behaving like a cosmological constant. This is shown by 
the subsequent oscillations of w with decreasing amplitude until w goes to —1. 

In the case of the Brane model the potential has a much sharper minimum. In that case the field stays longer in 
the attractor regime until suddenly it gets trapped in the minimum without oscillating. The effect is that w remains 
zero during most of the period of structure formation (with the field in the matter attractor regime) and then goes to 



1 This model was developed jointly with J. Weller [H 
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FIG. 2. Solutions obtained by including a V p factor in the potential 
radiation (a) (dotted). Radiation-to-matter transition is around a — 10 
Q.4, = 0.7 and Q, m = 0.247 today. 

LEFT: The AS model with A = 5, B = 54.4057, A = .01. Notice the behavior of the field at late times where fi 
near to zero. That is the period of structure formation so the field doesn't affect structure formation much. 
RIGHT: The Brane model with A = 5, B = 56.10425, A = .01, C = 1 and D = 0.1. Notice that the behavior of the field 
at late times is different from the AS model in the sense that the field retains a significant amount of energy density. This 
behavior affects structure formation as well as the CMB. 



0.0001 0.001 

a/a 

. We show f2^(a) (dashed), Q m atter(a) (solid) and 
~ 4 and a — 1 today. Both models have h = 0.65, 



goes very 



— 1 directly to start accelerating the universe. This allows the field to retain the significant amount of energy density 
seen previously. Moreover because of this, the Brane model starts accelerating the universe later than the AS model. 

Another question concerning both potentials above is the stability of the local minimum under quantum tunneling. 
It has been shown however that tunneling is negligible which renders the minimum in both potentials effectively 
stable g§. 

Finally we should say that the tightest constraint comes from requiring that not be too large during nucleosyn- 
thesis [§7| (at a « 1CT 10 ). For fi^(lAfeV) < 0.1 we get that A > 6.3. The models shown in Fig. | do not obey this 
constraint because they were chosen to dramatize the different behaviors. When comparing with data, we use more 
realistic parameters. 

III. DENSITY PERTURBATIONS 

To make connection with the real universe one has to consider the growth of small perturbations about an FRW 
metric. For extensive reviews see for example p^-|5l|. Here we adopt the conventions of Ma and Bertschinger fl5l| . 
We start with a line element given by 



a 2 (r)[-dr 2 + (<% + h i j)dx i dx j ] 



(13) 



The metric perturbation is split into scalar, vector and tensor modes as usual. 



A. Perturbations in the scalar field 



Here we give the relevant equations for the perturbations in the quintessence. For an extensive review of scalar 
field perturbations in cosmology see for example [pApOf and for the first applications in the case of quintessence P,p| . 
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FIG. 3. The equation of state parameter w (solid) and deceleration parameter q (dotted) for the same two models of Fig. ^. 
On the left we have the AS model and on the right the Brane model. The difference of the two models seen in Fig. ^ is seen 
again here in u>.(see the text for more details) Acceleration begins around a — 0.7 — 0.9 when q becomes negative. 



To include scalar fields into the perturbations lets write the field as <fr(x) = 4>o{ T ) + 4>i(s) with (f>i << <f>Q, where 
4>o(t) represents the zeroth order homogeneous field i.e. V0o = 0. We can then form the stress-energy tensor from 
which we can read the perturbed density pi, pressure Pi and velocity divergence 9 of the field. It is important to note 
that the shear is zero independently of the form of the potential. The quantities of interest in k — space are 



1 



Pi = -n<?Wl - V 4 



(14) 
(15) 
(16) 



The field perturbation also obeys a Klein-Gordon equation 



a ■ 
2-0i 
a 



+ (fc 2 + a 2 V M )(j)i + -h(j>o = 



(17) 



The above 2nd order equation is equivalent to the two 1st order equations (29) of [ pl| . One has to keep in mind 
however that in general(as noted also in |3(J), the density perturbations in the field are not adiabatic i.e. ^ is not 

in general equal to We shall use the symbols c w and c s for = ^ and c 2 s = the later being the adiabatic 
speed of sound. If the field is adiabatic then the two quantities become equal c w = c s . The above point is important 
for understanding why even if quintessence has the same equation of state as CDM it nevertheless does not clump. It 
all boils down to how the Jeans length is defined. In the case of CDM, the Jeans length is zero by definition, contrary 
to quintessence where it can be as large as the horizon. This is because the Jeans length must be defined in terms 
of c w and not c s to include the effect of entropy perturbations. In the case of quintessence it is exactly this effect 
which prevents clumping. Basically, the quintessence field is able to have a non-zero pressure perturbations even when 
c s = 0, and the pressure will resist gravitational collapse. 

Another thing to note about the field equation (JTt[) is the physical meaning of the various terms. The Hubble 
expansion provides a time-dependent "drag" term and the potential a time dependent "mass" term. The last term 
gives the gravitational sources for the field perturbations. 
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B. Growing super-horizon perturbations in the radiation era 



In our model, deep in the radiation era we can assume that the field is in the attractor and that the potential is 
a pure exponential. In this case we can use the analysis of Ferreira and Joyce |39| to get the initial conditions. We 
should note that the initial conditions given here are curvature initial conditions. 

We can use (10) with n — 4 for the radiation era. Then we have for Quintessence 



p<t> 


= w = -J 2 

4 


(18) 


a 
a 


1 

T 


(19) 




4 

~ A2 


(20) 




= 7j(t0i _ 01 ) 


(21) 




T(j>i + 01 
T(f>i - 01 


(22) 


T 2 (pl + 2t0i 


+ [(fcr) 2 +4]0j + = 
A 


(23) 



Assuming that only photons and quintessence are dominant during that period we find for the gravitational per- 
turbation h 

T 2 h + rh + 6Q 7 5 7 + 3 (! + 3c**)*V* = (24) 
For superhorison perturbations we assume kr << 1 and therefore to first order f? 7 = 0. Therefore 

Combining the last six equations (for k 2 = 0) we get 

T 2 h + Th-A{\ - ^)h+ |(2r0i +00 = (25) 
\ z A 

r 2 0i + 2t0! + 40i + - T h = (26) 
A 

Let h = Cr m and 0i = Dr m . Then 



-1±J§ -15 

m = ±2, m = * (27) 

The growing physical mode is the one proportional to t 2 . For the coefficients C and D we get (a result also found 
by Ferreira and Joyce) 

D = ~i c ^ 

This gives 

6* = J£ S c (29) 
Note that the above result gives = 3 7^ c 2 = w = i so quintessence is not adiabatic even if w = 0. 
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C. Initial Conditions 



From the above analysis we can form the initial conditions following the procedure in Ma and Bertschinger. The 
initial conditions for all the standard quantities remain the same. The two additional initial conditions for quintessence 
are 

1= -Afc, <t>x=-^h (30) 

Even if we don't use the correct initial conditions for the field, as a property of the growing curvature mode, the 
correct evolution of the field density contrast will be reached quite rapidly and the initial conditions given above will 
be valid. 

IV. STRUCTURE FORMATION 

First we investigate the growth of structure with time using the density perturbations as a guide. Then based on 
that investigation we give the resulting matter power spectrum and explain its form. 

A. The growth of structure 

One intuitive and important quantity is the dimensionless growth rate for the cold dark matter, 

n cS = T^. (31) 

This quantity give an instantaneous measure of the growth rate of structure. The time evolution of n e ff is shown for 
a number of models on the left panel of Fig. |]. Initially, outside the horizon, n e g = 2 completely independent of the 
model of structure formation(for a flat universe). This is due to the fact that because of causality the form of the 
background energy density is not relevant, only the amount. After entering the horizon, the growth of structure is 
suppressed during the radiation era simply because dark matter is extremely sub-dominant. During the matter era, 
dark matter becomes dominant again and in the case of an SCDM model, n c g eventually reaches the value of 2 again. 
For a ACDM model however, n e a after growing for a while, it eventually drops down to zero. This happens when the 
universe enters the accelerating era where the perturbations leave the horizon and structure stops forming. 

The two quintessence models considered here have a similar effect to a ACDM model. The only difference is the 
existence of a significant amount of dark energy during the matter era. This results in a suppression of structure 
formation since the amount of dark matter is less and therefore structure forms less efficiently. As stressed in section 
III. A, even when — quintessence can have pressure perturbations which allow it to resist clumping. 

This is not very important for the AS potential since if you notice from Fig. ^| £1^ goes very close to zero before 
the accelerating era. Therefore in that model structure is just a little bit suppressed compared to a ACDM model, 
but not by much. 

For the Brane potential however the above effect is much stronger. In fact stays quite significant during all of 
the matter era. Therefore for that potential, structure growth is even more suppressed as one can see from Fig. ^|. 
Not only do we have the acceleration era like before which suppresses structure anyway, we also have a significant 
amount of quintessence which doesn't cluster. This means that the amount of clustering matter is even less in this 
case so structure is even more suppressed. 

B. The matter power spectrum 

Based on the reasoning of the previous section, we can predict the form of the matter power spectrum. For the 
AS model we expect the power to be a bit less than for A model with the same cosmological parameters. The Brane 
model should exhibit additional suppression. This is indeed the case as one can see in the right panel of Fig. |[ 

If we use a Hubble constant of 75kms~ 1 Mpc~ 1 , even though structure is suppressed for the Brane model, it fits 
the data very nicely. All of the models shown in the figure have h = 0.75, fl c = 0.297, = 0.053, = 0.65 and 
some reionization with T opt = 0.1. The corresponding CMB anisotropies are shown later in the right panel of Fig. [t]. 
The C7 8 's 1.08,0.98 and 0.77 for ACDM, AS and Brane potentials respectively. 
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FIG. 4. LEFT: The effective dimensionless growth rate n e ff for four different models at k = 0.1Mpc~ . At early times 
n c ff — » 2 for all models, due to the r 2 dependence of the growing superhorizon modes. Growth is suppressed in the radiation 
era (t ps lOOMpc) and in the 0-era (today). In the matter era sCDM(dash-dotted) has the most growth, then ACDM (solid), 
AS (dotted) and finally with the least growth the Brane (dashed). The last two are the same models shown in Fig. ^. 
RIGHT: The matter power spectrum normalized to COBE using the Bunn-White fitting formula. The plotted models are A 
(solid), AS (dotted) Brane (dashed). All the models have h = 0.75, Q c ~ 0.297, fib = 0.053 and fi^, = 0.65. Here we use more 
realistic quintessence parameters consistent nucleosynthesis. In particular for the AS model A = 8, B — 33.9627, A = 0.01 and 
for the Brane model A = 8, B = 35.13689, C = 0.01 and D = 0.01. The data points are the decorrelated data of Hamilton 
et.al. [Eli - " 



It is not the goal of this paper to scan all of parameter space and compare a wide range of models, and we are not 
suggesting that any of the models in Fig. ^ are the best fitting model of their type. Figure |^ is intended to illustrate 
the physical differences of the different types of models by holding as many aspects as possible constant across the 
four models. Figure [| also illustrates the somewhat surprising result that it is possible for models with significant 
contributions from the dark energy at early times (i.e. the Brane model) to provide a good fit to the data. 



V. THE COSMIC MICROWAVE BACKGROUND 



Lastly we give the CMB anisotropy power spectrum and investigate the various physical effects due to the field. 
We follow the analysis of Hu and Sugiyama [ 
into primary and secondary. The primary anisotropies are the ones formed at the Last Scattering Surface (LSS), the 
secondary being the ones due to the subsequent cosmological evolution. Due to significant quintessence energy density 
at the LSS, the primary anisotropies are different compared with other models. In this section we drop sCDM from 
consideration and compare the remaining three models discussed above: the AS and Brane models, and the ACDM 
model with similar parameters. Moreover, since we are only interested in the differences cause by quintessence, we 
keep other cosmological parameters(f2;,, Q c and Hq) fixed. 



53- 56 1 . Following the common practice we divide the CMB anisotropies 



A. Primary anisotropies 



With the above parameters fixed, the only other parameter that can affect the primary anisotropies is the Hubble 
parameter during the time of LSS and earlier. Changing the Hubble parameter affects the heights of the peaks through 
three main effects, the driving effect, the doppler shift and Silk damping. 
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FIG. 5. LEFT: The Newtonian potential ip( T ) at fc = O.OlAfpc -1 , normalized to —1 as kr — > 0. Shown are the same three 
models of Fig. ^: AS model (dotted), A (solid), Brane model (dashed). Notice the very strong decay of the potential in the 
case of the Brane model, which leads to a very strong ISW effect. 

RIGHT: An illustration of the ISW effect for the same models. All models have the same initial power spectrum. The strong 
ISW for the Brane model at COBE scales affects the normalization and reduces the final anisotropies. 



Looking at fixed a it is easy to show that the Hubble parameter in a quintessence model Ha,{cl) is related to the 
Hubble parameter in a A model H\(a) at the same a, by 



Therefore, since for the quintessence models under consideration here is quite significant during all the history of 
the universe (where as Q\ isn't), the Hubble parameter in a quintessence universe will in general be larger than the 
one in a A universe. The above statement is of course true only if we keep flom, ^04> and Hq fixed which is what we 
do for comparing the models. 

The driving effect arises from a decaying Newtonian potential \& before the LSS (see Fig. as well as a dilation 
effect from a decaying space curvature $(not to be confused with <fr which is the quintessence field). The Newtonian 
potential and the curvature can be written in terms of the other perturbations as 

3a 2 H 2 v-^ r „ „ 3aH , 



2k 2 ^ L 1 k 2 

i 

&a^H 2 
k 2 



* = 3> ^— (fi 7 <7 7 + 3n„av) (33) 



The effect of quintessence on the driving effect comes from the significant fractional densities of quintessence during 
the radiation and matter era. During the radiation era fi^ is always greater than the on during matter era. This 
results in a greater difference of Hj, from H\ in the radiation era than in the matter era. The result is both the 
Newtonian potential and curvature decay faster in the case of quintessence. The effect is an increased driving effect 
in the case of quintessence which increases the temperature anisotropies for all peaks. We should also note that the 
same effect occurs by having more species of relativistic particles around(e.g. more massless neutrinos), the reason 
being exactly the same. Moreover this is more significant for the Brane potential whose Hubble parameter is always 
greater that the one for the AS potential. 

The velocity perturbation makes a contribution to the temperature anisotropies that is 90 degrees out of phase 
with the contribution from the density perturbation. The velocity is sub-dominant compared to the density which 
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FIG. 6. LEFT: The effect of Silk damping on the total primary anisotropics for the same models as Fig.|2| with the same 
initial power spectrum. In the case of quintessence, Silk damping is slightly more effective(see text). RIGHT: Primary 
temperature anisotropics for the same models, as would have been observed today if there was no other anisotropy source. The 
driving effect is observable on the first four peaks. The much higher 2nd peak for quintessence as well as the trough before it 
shows the Doppler effect. 



makes the effect weak. There is however a small but observable difference between a A model and quintessence. For 
quintessence even peaks as well as the troughs before them are slightly raised compared to a A model. This causes 
the oscillations appear to be "fattened" up due to the modulation of the velocity with the density to get the final 
temperature anisotropy. 

Finally we have Silk damping which is different in the three models. The damping coefficient ko is given by 



1 f a 1 R 2 + I(l + R) 
' L ' " !> ~6j a i H 2 r l (1 + R) 



kn(o,) = - „ 4Er2 _, n f m2 ' da (34) 



where R = and r opt is the optical depth. The damping coefficient is proportional to the Hubble parameter and 

is therefore smaller in the case of quintessence. This means that damping becomes more effective at smaller k (larger 
scales) in the quintessence models. The effect will be to suppress the small scale temperature anisotropies more 
effectively. Since at large k the driving effect eventually becomes unimportant, Silk damping eventually becomes 
stronger than the driving effect at small scales. This makes a A model with equivalent cosmological parameters have 
larger anisotropies than quintessence even though quintessence has larger anisotropies in the first two peaks. The 
total primary anisotropies are shown in Fig.|^ 

Let us now consider the secondary anisotropies which can change this picture quite a lot. 



B. Secondary anisotropies 

The most notable secondary anisotropy is the Integrated Sachs- Wolfe (ISW) effect. To dramatize this effect for 
pedagogical purposes we first consider the models from Fig. |^ which have a higher proportion of Quintessence than 
is allowed by nucleosynthesis (and an equally bad A model). 

The ISW comes from a decaying Newtonian potential after the LSS (see Fig. |5|). The standard treatment breaks 
the ISW down to early and late ISW. In our case however we have an intermediate ISW as well. 

The right panel of Fig. |jj shows ISW contribution for the three models while the left panel of Fig. |?] shows the 
angular power spectrum for the same models. The expression for Q lsw is the ISW contribution to the total anisotropies 
as defined in Hu and Sugiyama p5[ | eq.ll. 
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FIG. 7. The CMB anisotropies for the models of Fig.g. Shown on the left are the anisotropies with the same initial power 
spectrum and on the right after COBE normalizing. Due to the strong ISW the anisotropies of the Brane model (dashed) are 
suppressed after normalization. 

The ISW effect is very important for the Brane model (shown in all the relevant figures with a dashed curve). If 
what we saw today was just primary anisotropies then this model should have had the highest first peak of all three 
models. What we see however is exactly the opposite: it has the lowest peak. This is completely due to a very strong 
ISW which is shown on the right panel of Fig. Q The ISW for the Brane model boosts power at large scales^ < 100) 
compared to the other models. Therefore it will have a direct effect on COBE normalization (or in fact any other 
kind of normalization which includes points at those scales)- the anisotropies have to be scaled down to fit COBE. 
What happens physically is that to get the observed power, one has to use a smaller amplitude for the initial power 
spectrum. This makes the the ratio of the first peak to the spectrum at COBE scales smaller in the case of the Brane 
model. Normalizing to COBE then suppresses the small scale anisotropies. So even though the model had the highest 
peak at the LSS, it now has the smallest. 

The second effect of the very strong ISW on the Brane model is to shift the first acoustic peak to larger scales. 
This is though of lesser magnitude than the same effect caused by the smaller angular diameter distance. 

Finally we also have gravitational lensing on the CMB photons. This is more significant on smaller scales. In the 
case of polarization the effect is smaller than 10% and for the temperature anisotropies even smaller | p7| . The effect is 
a smearing of the peaks due to mixing of Z-values without destroying the overall structure of the peaks. We calculate 
it using the linear evolution method of CMBFAST and assuming that quintessence will have the same effect on lensing 
as a pure cosmological constant. This is reasonable since like a cosmological constant, quintessence doesn't cluster 
and therefore its power spectrum does not contribute to lensing. The quintessence power is actually much smaller 
than the photon or neutrino power which themselves make a vary small contribution and are neglected. 



The position of the peaks depends on two quantities, the sound horizon and the angular diameter distance to LSS. 
Both are affected by the presence of quintessence. 

First lets consider the sound horizon. Even though the baryon density is kept the same, hence the speed of sound 
is the same in all three models, the time to the LSS is different due to the different expansion rate. Hence the sound 
horizon given by 



C. Position of the peaks 




(35) 
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FIG. 8. LEFT: Tensor C;s for cosmological constant and the two quintessence models. All models have h — 0.65, fit, = 0.053 
fi c = 0.247, Q$ — 0.7, n s — 0.98, n t = —0.02 and ratio of tensor to scalar quadrupole of 0.1. Shown is a ACDM(solid), AS 
(dotted) with A = 5, B = 54.4057, A = .01 and Brane (dashed) model with A = 5, B = 56.10425, A = .01, C = 1 and D = 0.1. 
RIGHT: The CMB anisotropics for the more realistic models of Fig . ^|(Right panel). The effects discussed in the text are still 
visible but not as strong. Shown are the BOOMERanG98(solid) (Hsj and MAXIMA(dashed) ^ data. 

is different (t* is the time at last scattering). More specifically it is inversely proportional to H which means that 
quintessence has a smaller sound horizon. This shifts the anisotropies at LSS(in fc-space) to smaller scales. Greater 
H also makes the time at LSS r* smaller. 

Next we have the angular diameter distance to LSS which is smaller in the case of quintessence. For the AS model 
it is not much smaller than a A model since the universe starts accelerating at around the same time in both models. 
In the case of the Brane model however that universe is younger and therefore the angular diameter distance to 
LSS, smaller. In fact the more quintessence we have during the history of the universe, the younger it is since the 
universe starts accelerating later. If we have enough quintessence it doesn't accelerate at all. This appears to be 
counterintuitive. Remember however that in order to have more quintessence around one has to make the potential 
shallower. If it is shallow enough the field just rolls through. 

From the analysis of ]55fl the difference in I between peaks is given by 

SI = Tr^Il (36) 
r s 

where To is the time today. For the models of the right panel of Fig. |^ we get 292, 297 and 293 which is in very good 
agreement with the actual numerical result. The AS model has roughly the same angular diameter distance as the 
ACDM model but a smaller sound horizon. This makes the peak separation larger. The Brane model has the smallest 
sound horizon of all, but it also has the smallest angular diameter distance to the LSS since the universe is younger. 
The two roughly cancel each other and give roughly the same peak separation as in a ACDM model. 

Then we have the position of the first acoustic peak. The AS model has a first peak roughly at the same scale 
as the A model. The Brane model however has the first peak shifted to larger scales. This is a combination of two 
effects. On is the much smaller angular diameter distance and the other is the very strong ISW. 

D. Tensor modes 

As mentioned previously, the scalar field has no shear and therefore the only effect it has on the tensor modes is 
through the background evolution. The tensor anisotropies are shown in FigM 
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FIG. 9. LEFT: E-type polarization for the same models as Fig. [| The strong ISW in the case of the Brane model suppresses 
indirectly the polarization spectrum. 

RIGHT: Tensor E-type(heavy curves) and B-type(light curves) polarization for the models same models as Fig.[s](left panel). 
Since there is no ISW the difference in the AS and A models is just in the frequency of oscillation. For the Brane model we do 
have however the indirect effect of a strong ISW through normalization of the scalar Cis amplitude 



1500 



The amplitude of the tensor modes is in general much less than the scalar amplitude and the only important 
contribution comes at large scales. Even though they don't affect the temperature anisotropy much they are still 
important as they create B-type polarization which might be observed. 

The main effect is the strong ISW of the Brane model which suppresses the anisotropies at very large scales through 
the mechanism discussed in the earlier subsection. In the case of the AS model, the ISW is not as strong, so the model 
still retains anisotropies of roughly the same amplitude as a A model. Finally we have the shifting of the oscillation 
pattern due to the different sound horizon and angular diameter distance just like the scalar case. 



E. Polarization 



Polarization depends only on the physics at the LSS, unlike the temperature anisotropy which also depends on the 
evolution since the LSS. It is therefore a direct probe of the physics at the LSS. For a review of polarization see |Q 

E-type polarization for scalar modes is shown in Fig. |9|. The first feature we observe is a difference in the position 
of the peaks. This is due to the same effects as in the scalar C\ case, i.e. the sound horizon and angular diameter 
distance. 

The heights of the peaks are affected by the driving effect on the local quadrupole at LSS. Therefore in the left panel 
of Fig. U we see again as expected that the AS model has the largest amplitude followed by the A model. The Brane 
model should have had the largest amplitude but we have again the indirect effect of the ISW on the normalization 
which suppresses the whole thing. This happens again in this case because we normalized to COBE- the mechanism 
of suppression being the same as the one described in the scalar Ci case. 

B-type polarization is shown on the right of Fig|]. To get a sizable B-type polarization we used slightly different 
cosmological parameters(shown in the caption) from the ones used so far, basically changing the spectral indices and 
having a non-zero tensor to scalar quadrupole amplitude. The difference in the AS and A models is just in the frequency 
of oscillation. For the Brane model we do have however the indirect effect of a strong ISW through normalization of 
the temperature anisotropies amplitude (same thing as in the E-polarization). We should note however that lensing 
can also produce B-type polarization but we don't discuss it here since in the case of these models quintessence can 
be treated like a cosmological constant as far as lensing is concerned. 

Finally we have the polarization-temperature cross correlation shown in Fig.[lO|. The importance of the cross 
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FIG. 10. Cross correlation of temperature and polarization. On the left we have scalar cross correlation for the same models 
as Fig. ^. On the right we have tensor cross correlation for the same models as Fig|| In both cases we see a strong suppression 
of the correlation amplitude for the Brane model due to the ISW effect. 

correlation has been discussed in |n| for scalar modes and |32] for tensor modes. Again we see the different oscillation 
frequencies due to the different sound horizons and angular diameter distances just like the temperature anisotropies 
case, for both the scalar and tensor modes. In the case of the Brane model the amplitude of the correlation is smaller 
simply because in this model the temperature anisotropies (after normalization) are smaller compare to the other 
models. Have we kept the same initial power spectrum for all models the Brane model would have had the highest 
correlation. In the scalar case, on scales between I > 200 and I < 1000 we have the driving effect again giving more 
correlation to the quintessence models. Moreover for I < 600 we also have the ISW boosting the correlation for 
quintessence. At very small scales the correlation amplitude tends to be the same in all models. 



Before concluding we comment briefly on the magnitude-redshift relation for the brane model. The magnitude- 
redshift relation m(z) for a standard candle of fixed fiducial absolute magnitude is an essential tool for pinning down 
the nature of the cosmic acceleration. Different quintessence models will produce different functions m(z) which can 
be discriminated using data from, for example, type la supernovae (see for example p3|]). When considering possible 
future probes of m(z) it is essential to subject the space of possible models to the constraint that the models give a 
realistic account of cosmic structure. In particular, one can expect the structure-suppressing nature of quintessence 
to result in substantial constraints on the significance of quintessence at z greater than around unity since structure 
needs to have a chance to form |64|]. It is therefor intriguing that that the Brane model (which gives a good fit to 
cosmic structure data) has some pretty interesting features in m(z). 

Figure [ll] shows Am(z) for the Brane model (identical to the model shown in figure ||), a model for which the 
dark energy has an equation of state p — wp with constant w=-0.9559, and simulated binned data representing 1900 
supernova events as might be provided by the proposes SNAP satellite |6^| . The m(z) curve for a pure cosmological 
constant model (which was used to produce the simulated data) is subtracted from m(z) for the other models to get 
Am(z). Except for the choice of theoretical models shown, Fig. [ll] is identical to Fig. 2 in Q, where details of its 
construction can be found. 

The Ato(z) curves for the two dark energy models are essentially identical up to z « 1.3 at which point they 
diverge. As discussed in |33],[35| there is no question that a SNAP-type data set will have a high impact on our ability 
to discriminate among models of the dark energy. Our work on the Brane model illustrates how even with a high 
quality SNAP-class data set in hand, there could be exciting opportunities to further discriminate among realistic 
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models if a standard candle could be found that is effective out to higher z. 




FIG. 11. Am(z) for the Brane model (heavy curve) and atc = const. = —0.9559 model (dashed curve), along with simulated 
SNAP data. The curves are essentially identical up to z ~ 1.3 at which point they diverge. If a standard candle could be found 
that is effective at high redshifts it could play a significant role in discriminating among these types of models (We thank J. 
Weller for producing this figure). 



VII. CONCLUSIONS 

We have presented a detailed investigation of a class of quintessence models motivated by our earlier work [ [L9[ . 
These models employ a particular mechanism (the exponential potential with a prefactor) which allows realistic 
models to be produced with all potential parameters O(l) in Planck units. The potentials have a reasonable chance of 
developing strong theoretical foundations in brane theory or other theories with extra dimensions. The work presented 
here has taught us that these models also have an interesting and potentially observable impact on the formation of 
cosmic structure. The same mechanism that makes these models attractive from the point of view of fundamental 
physics causes the quintessence to play a much more significant role throughout the history of the Universe. This 
feature leads to interesting effects on the microwave background anisotropics, the matter power spectrum and the 
magnitude redshift relation that result in potentially observable differences from the predictions of other dark energy 
models. We have extensively examined the physical causes and the nature of these effects. 



To illustrate the possibilities, Fig. 12 shows the full CMB Temperature anisotropy power spectrum for the Brane 
(right) and AS (left) models, containing all the effects discussed in this paper. Each plot also shows the A model which 
best mimics the corresponding quintessence model. The lower panel of each plot shows that despite the closeness of 
the two curves, they are potentially differentiable by the MAP data (and of course we will eventually have much more 
than just the MAP data set). 

One issue which still needs to be addressed is the degree to which the signals we have discovered can realistically be 
differentiated from all possible mimicking behavior due to the dependence of the predictions on a range of cosmological 
parameters (the mimicking A models in Fig. [l^ were produced "by hand" , by making a thorough but not completely 
exhaustive exploration of all possible parameters). This paper lays the groundwork for such a project. The full impact 
of our results will not be know until this this issue is addressed in a more systematic way. 

Still, it is quite interesting that our quintessence models give predictions that fit all existing constraints, and which 
leave a potentially unique set of signals that could be observed with future experiments. 



1G 




500 1000 
multipole moment, I 



1500 



4 




500 1000 
multipole moment, I 



FIG. 12. A possible future comparison of the AS and Brane models with a A model using MAP data. On the left we have 
the AS model of the right of Fig.§ compared with a A model with h = 0.75, tt A = 0.7005, tt b = 0.053, tt c = 0.2465 and r = 0.1. 
On the right we have the Brane model of the right of Fig.^ compared with a A model with h = 0.75, J7a = 0.66327, Sit = 0.053, 
fl c = 0.28373 and r op t = 0.1. The two A models were chosen to match the first peak and at the same time the rest of the curve 
as close as possible. On the bottom of each figure we show the residuals of the two corresponding models with MAP error-bars 
(We thank L.Knox for providing the MAP error-bars). 
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